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Deleptonization and a Nuclear Equation of State 
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We present results from the first 2D and 3D simulations of the collapse of rotating stellar iron 
cores in general relativity employing a finite-temperature equation of state and an approximate 
treatment of deleptonization during collapse. We compare fully nonlinear and conformally flat 
spacetime evolution methods and find that the conformally flat treatment is sufficiently accurate 
for the core-collapse supernova problem. We focus on the gravitational wave (GW) emission from 
rotating collapse, core bounce, and early postbounce phases. Our results indicate that the GW 
signature of these phases is much more generic than previously estimated. In addition, we track the 
growth of a nonaxisymmetric instability of dominant m = 1 character in one of our models that 
leads to prolonged narrow-band GW emission at ~ 930 Hz over several tens of milliseconds. 
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Introduction. — For more than two decades astrophysi- 
cists have struggled to compute the gravitational wave 
(GW) signal produced by rotating stellar core collapse 
and the subsequent supernova evolution. Besides the co- 
alescence of black hole and neutron star binaries, core- 
collapse events are considered to be among the most 
promising sources of detectable GWs. Theoretical pre- 
dictions are still hampered by three major problems: (i) 
The rotational configuration prior to gravitational col- 
lapse is still uncertain since multi-D evolutionary calcu- 
lations of rotating massive stars have not yet been per- 
formed; (ii) reliable waveform estimates require a general 
relativistic (GR) treatment, since both high densities and 
velocities in combination with strong gravitational fields 
are encountered in this problem; and (iii) an adequate 
treatment of the nuclear equation of state (EOS) and 
the neutrino microphysics/radiative transfer is crucial for 
obtaining realistic collapse, bounce, and postbounce dy- 
namics and waveforms. GW emission from core-collapse 
supernovae may arise from rotating collapse and bounce, 
postbounce neutrino-driven convection, anisotropic neu- 
trino emission, nonaxisymmetric rotational instabilities 
of the protoneutron star (PNS), or from the recently pro- 
posed PNS core g-mode oscillations. Previous estimates 
of the GW signature of core-collapse supernovae have ei- 
ther relied on Newtonian simulations [l ,J1, H, 0, H, [B| (to 
some extent approximating GR effects |7[) or GR simu- 
lations with simplified analytic (so-called hybrid) EOSs 
and no neutrino treatment 

BBEIHl- Depending on 
the rotation strength, softness of the EOS at subnuclear 
densities, and inclusion of GR effects, the collapse dy- 
namics and accordingly the GW signature can differ sig- 
nificantly. 

Here we present new results from GR simulations, fo- 
cussing on the rotating collapse, bounce, and early post- 
bounce phases. These are the first-ever multi-D simula- 



tions in GR with presupernova models from stellar evolu- 
tion calculations, a finite-temperature nuclear EOS, and 
a simple, but effective treatment of electron capture and 
neutrino radiation effects during collapse. In this way we 
obtain the most accurate estimates of the GW signature 
of rotating stellar core collapse in full GR to date. 

Method and Initial Models. — We perform all 3D sim- 
ulations in full 3-1-1 GR using the Cactus infrastruc- 
ture Cartesian coordinates, and mesh refinement 
provided by the Carpet driver [l3|. The only assump- 
tion on symmetry is refiection invariance with respect 
to the equatorial plane. Spacetime is evolved using the 
BSSN formulation (see, e.g., [l3|) and we fix the gauge 
freedom by 1 -f log slicing and by a hyperbolic shift[lf!|. 
We use the hydrodynamics code Whisky QJj] , which im- 
plements the equations of GR hydrodynamics via finite- 
volume methods. Typical simulation grids extend to 
3000 km and use 9 refinement levels. The central res- 
olution is 350 m. In addition, we perform axisymmet- 
ric (2D) simulations for all models using the CoCoNuT 
code [a, [l3| , which approximates GR by the conformal 
flatness condition (CFG). CoCoNuT utihzes spherical 
coordinates with 250 logarithmically spaced radial and 
45 equidistant angular zones, covering 90°, and a central 
resolution of ^ 250 m. Resolution test calculations with 
both CoCoNuT and CCW do not yield significant quali- 
tative or quantitative differences. We extract GWs using 
a variant of the Newtonian quadrupole formula 91 . 

We employ the finite-temperature nuclear EOS of Shen 
et al. Il8l| . Deleptonization is implemented as proposed 
by [l9|: The electron fraction Y^, is parameterized as a 
function of density based on data from spherically sym- 
metric radiation-hydrodynamics calculations with stan- 
dard electron capture rates '^2^. After core bounce, Yg 
is passively advected and further lepton loss is neglected, 
but neutrino pressure contributions continue to be taken 
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TABLE L Summary, pb is the density at bounce, the max- 
imum characteristic GW strain ftchar.max is at a distance of 
10 kpc, and iSgw is the energy emitted in GWs. Values for 
E20Apb include GW emission from late-time 3D dynamics. 
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into account above trapping density [l9|. 

In this study we focus on the collapse of massive pre- 
supernova iron cores with at most moderate differen- 
tial rotation, and rotation rates that may be too fast to 
match garden- variety pulsar birth spin estimates [2l], , 
but could be relevant in the coUapsar-type gamma-ray 
burst context. As initial data we use the non-rota ting 
20 Mq presupernova model s20 of Woosley et al. [23| 
which we force to rotate according to the rotation law 
discussed in [1, 0] . We parameterize our models in terms 
of the differential rotation parameter A and the initial 
ratio of rotational kinetic to gravitational binding energy 
j3i = T /\W\. In addition, we perform a calculation with 



the 20 Mq model E20A of Heger et al. [2J], which was 
evolved to the onset of collapse with a ID treatment of 
rotation. In Table[T]we summarize the model parameters. 

Results. — In Fig. [1] we compare GW signals computed 
with CoCoNuT in 2D-CFC and those computed with 
our 3D-full-GR approach. Model s20A2B2 (red hues) is 
a moderate rotator with a /3i = 0.50%, rotating rigidly 
in its central region. It stays axisymmetric throughout 
its numerical evolution. The agreement of 2D-CFC with 
3D-full-GR is excellent for this model: Both waveforms 
match almost perfectly at bounce and during the very 
early postbounce phase. A few ms after bounce, when 
convection in the region behind the stalling shock sets 
in due to a negative entropy gradient, the signals begin 
to differ quantitatively while remaining in phase. We 
attribute this small mismatch to the choice of coordinate 
grids and to differences in the growth and scale of vortical 
postbounce motions between 2D and 3D. Model s20AlB5 
rotates with constant in the entire core. Despite its 
very large (3i = 4% it remains essentially axisymmetric 
during the time covered by our simulation, since most of 
its angular momentum is attached to material at large 
radii that falls inward and spins up only slowly. The 
waveforms in CFC and full GR agree very well. Again, 
both waveforms match best for the strong burst related to 
core bounce during which more than ~ 90% of the total 
GW energy are emitted in an axisymmetric model. The 
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FIG. 1: GW strain h+ along the equator for models s20A2B2 
and S20A1B5. We compare 2D-CFC and 3D-full-GR results. 



6 

4 

. 2 

-2 
-4 
, -6 
' -8 
-10 
-12 





1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 . 










S20A2B1 " 






S20A2B2 -: 






S20A2B4 J 






S20A3B3 J 






E20A ; 







■ 



t-tb 



, (ms) 



10 



15 



FIG. 2: Equatorial GW strain h+ for representative models. 



overall excellent agreement of CFC with full GR confirms 
results of [1, [HI and proves that CFC is a very good 
approximation to full GR in the core-collapse scenario. 

In Fig. [2] we present waveforms of models with vary- 
ing initial degree of differential rotation A and rotation 
rate /3i. We find that the inclusion of a microphysical 
finite-temperature EOS and of electron capture during 
collapse yields results that differ considerably from those 
obtained in previous, less sophisticated studies. Fig. [2] 
exemplifies that largely independent of the initial rota- 
tional configuration, the GW signal of rotating collapse 
has a generic shape: a slow rise in the pre-bounce phase, 
a large negative amplitude at bounce when the motion 
of the inner core is reversed, followed by a ring-down. 
This so-called "Type I" signature corresponds to a bary- 
onic pressure dominated bounce [l|, [3, 0,1^ • Thus all our 
models undergo core bounce dominated by the stiffening 
of the EOS at nuclear density. 

This is in stark contrast to the studies using the hy- 
brid EOS 0, i, i, Ei, where initial models with rota- 
tion rates in the range investigated here develop sufficient 
centrifugal support during contraction to stop collapse 
at subnuclear densities, resulting in several consecutive 
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FIG. 3: Normalized mode amplitudes Am at postbounce 
times (upper panel) and GW strains h+ and hx along the 
poles (lower panel) for model E20A. 



centrifugal bounces separated by phases of coherent re- 
expansion of the inner core. While in GR models exhibit- 
ing such multiple centrifugal bounce and the correspond- 
ing "Type 11" GW signals are only rarer compared to 
Newtonian gravity Q , we do not observe any such model 
in this study. An evident example is model s20A2B4: In 
previous studies without detailed microphysics, the cor- 
responding model with identical initial rotation parame- 
ters (A2B4G1) showed clear "Type 11" behavior in both 
Newtonian and GR calculations 0, Q . 

The suppression of the multiple centrifugal bounce sce- 
nario is due to two physical effects: (i) In contrast to the 
simple hybrid EOS, in our case the mass and dynamics 
of the inner core (which is most important for the GW 
emission) is not merely determined by the adiabatic in- 
dex 7 = d\nP/d\np (at constant entropy) of the EOS, 
but also by deleptonization during collapse. This leads 
to considerably smaller inner cores with less angular mo- 
mentum and weaker pressure support, (ii) Since multiple 
centrifugal bounce was observed for a model with ini- 
tially moderately fast rotation in a previous Newtonian 
study [l[ where both a microphysical finite-temperature 
EOS and a deleptonization scheme were employed, the 
absence of this collapse type in our study is not only 
caused by microphysical effects, but also by the effec- 
tively stronger gravity in GR. This is in accordance with 
simulations using the simple hybrid EOS [8]. A detailed 
analysis of the interplay and quantitative influence of 
the above two effects responsible for the elimination of 
multiple centrifugal bounce in the rotational stellar core- 
collapse scenario will be presented in a future publication. 

Model E20A possesses the largest GW amplitude of 
all our models. In addition, it reaches a high /3b at core 



bounce (see Table IJ) and settles at a postbounce (3{ of 
^ 9%. In a previous Newtonian study, [25[ have found a 
low-T/jM^I nonaxisymmetric instability for a PNS with 
similar /3f. In order to verify their findings, we trace the 
evolution of model E20A to 70 ms after bounce and per- 
form an analysis of azimuthal density modes oc e™*^ in 
the equatorial plane by computing complex Fourier am- 
plitudes Cm — 57 J^^ p{w,ip, z = 0)6""*^ dtp on rings of 
constant coordinate radius with respect to the coordi- 
nate center of mass. The latter stays within the inner- 
most zone at all times. In the top panel of Fig. [3] we 
display the normalized mode amplitudes Am = ^'%~ 
tracted at 15 km radius. Without adding artificial seed 
perturbations to model E20A, discretization errors trig- 
ger m = {1, 2, 3} modes, which rise to a level of ~ 10~^ 
during the ^ 220 ms collapse phase. After bounce, the 
m=l mode exhibits the fastest growth. This growth on 
a dynamical time scale, lasting over tens of ms until sat- 
uration, is closely followed by a growth of m — {2, 3} 
daughter modes 2^ 2^ . Note that the m — A Cartesian 
grid mode starts out on the ~ 10"'* level and remains 
constant until the plunge phase of collapse during which 
all modes are amplified. After core bounce, model E20A 
remains dynamically stable to the m = 4 grid mode. In 
the lower panel of Fig. [3] we plot the GW strains /i+ and 
hx as seen along the polar axis. The rotational symmetry 
of E20A before and early after bounce times is reflected 
in the fact that h+ and hx are practically zero until E20A 
develops considerable nonaxisymmetry at ~ 25 ms after 
bounce when the m = I mode becomes dominant and its 
TO = 2, GW-emitting harmonic reaches a sizable ampli- 
tude. In remarkable agreement with expectations for a 
spinning bar, and hx oscillate at the same frequency 
(~ 930 Hz) and are phase-shifted by a quarter cycle. 

Discussion. — Our results indicate that the GW signa- 
ture of the collapse, bounce, and early postbounce phases 
of the core-collapse supernova evolution is much more 
generic than previously thought. We find that the dy- 
namics of core bounce are dominated by mainly gravity 
and microphysics, reducing the importance of centrifugal 
support for the rotation rates considered here. Impor- 
tantly, for our model set we do not observe rotationally 
induced multiple core bounce at subnuclear density as 
proposed by previous studies that did not include a mi- 
crophysical finite-temperature EOS and electron capture 
treatment in combination with GR. Thus we predict that 
the core-bounce waveform of models in a large parame- 
ter space of initial rotation rate and degree of differential 
rotation will likely both qualitatively and quantitatively 
resemble those presented in Fig. [5] 

Model E20A, which we evolve to later postbounce 
times, exhibits the dynamical growth of a nonaxisym- 
metric low-T/|iy| corotation-type to = 1 instability 
H, [ill. We also find m = {2, 3} contributions and 
thus significant GW emission from the quadrupole com- 



4 




E. Seidel, D. Shoemaker, N. Stergioulas, and J. Tohline 
for help and stimulating discussions. This research was 
partially supported by the DFG (SFB/Transregio 7) and 
by the NCSA under grant No. AST050022N. 



f(Hz) 



FIG. 4: Spectra of the characteristic GW strain /ichar of all 
models and the LIGO (optimal) rms noise curves [30| . 



poncnts of these modes. We emphasize that we observe 
this instability not only in E20A, but also in other mod- 
els with comparable values of /3f. Our results, which re- 

0,[25, 



28, 2 



move the limitations of previous studies 
demonstrate that the development of nonaxisymmetric 
structures is neither limited to Newtonian gravity, sim- 
ple matter models, equilibrium configurations, nor high 
values of /3, but may rather be a phenomenon occurring 
generically in differentially rotating compact stars. 

For the GW signals from the axisymmetric collapse 
and core-bounce phase we obtain peak amplitudes of up 
to h ^ 10"^*^ at 10 kpc, while the nonaxisymmetric struc- 
tures in model E20A developing later emit GWs with 
h only a factor ~ 5 smaller. However, since the latter 
emission process operates over several tens of ms, the 
total energy Eg^ emitted in GWs is larger than that 
from the core-bounce signal. This is evident in Fig. 31 
where we display the characteristic GW strain spectra 
^char = ^y2^^^'^dEgv,/df [5] for all models, evolving 
E20A for 70 ms after bounce (see also Table IJ). Consid- 
ering only the core-bounce waveforms, /ichar has its max- 
imum between 300 and 800 Hz, while for model E20A it 
peaks at ~ 930 Hz, the frequency of the GW-emitting 
component of its nonaxisymmetric structures. We con- 
clude that the core-bounce GW signals of all models in- 
vestigated here may be detectable by current and future 
LIGO-class detectors from anywhere in the Milky Way. 
Models that develop nonaxisymmetric instabilities may 
be detectable out to much larger distances if the insta- 
bility persists for a sufficiently long time. 

We point out that due to the nature of the approxi- 
mation used for the neutrino effects in this study, we can 
only accurately model the GW emission in the collapse, 
bounce, and early postbounce epoch of the core-collapse 
supernova scenario, but not much later than the neu- 
trino burst at shock breakout a few ms after bounce. In 
future work we plan to improve upon this and carry out 
longer-term postbounce evolutions, where additional GW 
emission mechanisms may play an important role 0] • 
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